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PROPELLANT VAPORIZATION AS A CRITERION FOR ROCKET-ENGINE DESIGN- 
CALCULATIONS OF CHAMBER LENGTH TO VAPORIZE VARIOUS PROPELLANTS 

By Richard <J. Priem 


SUMMARY 

Vaporization rates were calculated for drops of n-heptane, ammonia, 
hydrazine, oxygen, and fluorine. The percent propellant vaporized is 
correlated with an effective chamber length for various spray conditions, 
and various engine-design and operating parameters. The results show 
that the effective chamber length required to vaporize a given high per- 
centage of propellant is the shortest with oxygen and increases for 
fluorine, heptane, ammonia, and hydrazine in that order. 


INTRODUCTION 

Calculations were reported in references 1 and 2 for the rate of 
fuel vaporization in the combustion chamber of a n-heptane - oxygen rocket 
engine. How the vaporization-rate calculations could be used to predict 
combustion efficiencies and to design combustors was also indicated. 

The results of these calculations, based on a combustion model in which 
vaporization of the fuel was rate controlling, showed how various design- 
and operating-parameter changes would affect the vaporization rates of 
n-heptane drops burning in oxygen. Reference 2 also indicated that a 
small number of large drops that do not vaporize completely may be re- 
sponsible for much of the loss in rocket -engine performance. Experimental 
results obtained with an engine (ref. 3) using n-heptane as the fuel 
agreed with the calculations for sprays having geometric standard devia- 
tions of 2.5 and mass median drop radii of 70 to 280 microns depending 
on the type of injector. 

This report covers the additional calculations made at the KACA 
Lewis laboratory to determine the vaporization rates of oxygen, fluorine, 
ammonia, and hydrazine in rocket engines. The results are presented in 
terms of an effective chamber length required to vaporize a given percent 
of the propellant . The analysis assumes that the vaporizing propellant 
is burned instantaneously at stoichiometric conditions. The analysis 
further assumes that there is one -dimensional steady-state flow and that 
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all drops are produced at the injector face and have equal velocities. 
The model and equations of reference 1 are used to determine the veloc- 
ity , mass, and temperature of the individual" drops at various distances 
downstream. 


METHOD 


The calculation used the model, the iteration technique, and the 
following equations from references 1 and 2: 


Mass -transfer rate, 


Heat -transfer rate, 


Drop-heating rate, 


Drop acceleration, 


Gas velocity. 
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For calculation purposes the propellant was assumed to he in five drop- 
size groups, which were arbitrarily chosen from equation (6) as those 
equal to 10, 30, 50, 70, and 90 percent of the mass in drops smaller 
than the drop radius. With a median drop radius of 75 microns and a 
standard deviation of 2.3 the five .drop radii selected by this method 
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were 25, 48, 75, 120, and 225 microns . The number of drops in each group 
was so chosen that each group had 20 percent of the total mass. The 
symbols used in this report are defined in appendix A. 

The equations were solved by an iterative procedure as described 
in reference 4. The average physical properties in the mantle surround- 
ing the drop were also determined by the method described in reference 4. 
Table I shows the design and operating conditions used in the calculations. 
Appendix B shows equations used to describe the physical properties for 
heptane, ammonia, and hydrazine vaporizing in a gaseous -oxygen atmos- 
phere, oxygen vaporizing in a gaseous -heptane atmosphere, and fluorine 
vaporizing in a gaseous -hydrogen atmosphere. 


RESULTS MD DISCUSSION 

The results of the calculations are droplet histories, examples of 
which are shown in figure 1. Droplet temperature, droplet velocity, 
gas velocity, percentage of mass vaporized, and vaporization rate are 
shown at various chamber lengths for n-heptane, ammonia, hydrazine, 
oxygen, and fluorine. The conditions used in the calculations shown 
in figure 1 are indicated in table I. The initial temperature of heptane 
and hydrazine was taken at room temperature. Ammonia was arbitrarily 
taken at its normal boiling point while the initial temperature for 
fluorine and oxygen corresponds to the normal boiling point of nitrogen. 
The effect of selecting other temperatures and conditions is indicated 
in reference 1 and later in the Correlation of Results section. 


Temperature Histories 

The temperature histories of the small drops (25 microns), median 
drops (75 microns), and large drops (225 microns) in the spray are 
shown in figure l(a) for five propellants. The calculated temperatures 
rise to a steady value corresponding to a wet-bulb temperature. The 
wet-bulb temperature and corresponding vapor pressure are tabulated in 
the following table for the various propellants. Also tabulated is the 
ratio of the distance required to attain the wet -bulb temperature to 
the distance required to vaporize the drop. 


Propellant 

Wet -bulb 
tempera- 
ture, 

°R 

Vapor pressure 
at wet bulb, 
Ib/sq in. 

Length to wet bulb 
Length to vaporize 

Heptane 

845 

133 

V . 6 

Hydrazine 

860 

165 

1/15 

Ammonia 

555 

205 

1/15 

Oxygen 

235 

275 

1/15 

Fluorine 

220 

255 

1/15 
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Gas -Velocity Histories 

The average velocity of the vaporized and "burned propellant at 
various positions down the chamber is shown in figure l(b) . For all 
propellants this gas velocity initially increases rapidly and then 
asymptotically approaches the final gas velocity. The gas -velocity 
curves for oxygen and fluorine are about the same (liquid oxygen being 
slightly higher than fluorine) and were 1 higher than the other three 
propellants. The lowest gas velocity curve was obtained with hydrazine. 
The ranking order of propellants with respect to decreasing gas veloc- 
ity at a fixed chamber length was (l) oxygen, (2) fluorine, (3) n-heptane, 
(4) ammonia, and (5) hydrazine. 


Droplet -Velocity Histories 

Droplet -velocity curves for the small-, median-, and large -diameter 
drops are shown in figure l(c) for the five propellants. The velocity 
of the drop when it is 99 percent vaporized^is indicated by the solid 
symbol. The lowest velocity at complete vaporization was obtained with 
the high-density fluorine. The small drops accelerate faster than the 
larger drops; thus, a mixing effect occurs because of the relative motion 
of the drops. 


Mass Histories 

Curves showing the percent mass vaporized of the small, median, and 
large drops in the spray as well as the average for the spray are shown 
in figure l(d) for the five different propellants. The order of the 
propellants, based on vaporization rate; was (l) oxygen, (2) fluorine, 
(3) n-heptane, (4) ammonia, and (5) hydrazine. The shapes of all the 
curves were similar with the exception of n-heptane which crosses the 
ammonia curves. The curves show that in 2 .inches the small drops of 
all propellants are almost completely vaporized. However, in 2 inches 
only a small percentage of the large drops "was vaporized. 


Vaporization-Rate Histories 

The total-vaporization-rate per unit-length curves (fig. 1(e)) are 
similar to those shown in references 1 and 2. There are two peak points 
and a minimum for each propellant. The highest peak vaporization rate 
is obtained with fluorine. Second highest peak is with oxygen, third 
with ammonia, fourth with heptane, and the lowest peak occurs with hydra- 
zine. The peak values occur at different positions in the chamber with 
the various propellants. The peak occurs closest to the injector face 
with oxygen and fluorine. The peak is the farthest from the injector 
with heptane and hydrazine. Ammonia falls between the two extremes. 


4720 



4/2U 


NACA TN 3883 


5 


Correlation of Results 

The calculated, results obtained for the conditions listed in table X 
were correlated for each, propellant . Percent mass vaporized and unvapor- 
ized curves for sprays having a standard deviation Cg of 2.3 were 
plotted against effective length in figure 2 for each propellant. This 
deviation represents a spray in which 68.23 percent of the mass is in 
the drop-size range between Mg yj a g and CgMg The correlation of 

results was accomplished by modifying chamber length with a factor to 
obtain an effective length. 


Effective length = 


LP° . 66 u 0^40 1 _ 9x]_Q "5 

7 „ vO. 40.. 1.45 0.75 

(1 - ng Mg^ M v Q 


This factor was obtained from crossplots of the lengths required to 
vaporize 90 percent of the mass. The correction factor is slightly dif- 
ferent from that determined in references 1 and 2. The critical tem- 
perature was used to obtain a reduced initial temperature T r . The term 

(l - T r ) was used to indicate that the model of vaporization does not 

apply when the initial temperature is greater than the critical tempera- 
ture (in this situation T r > 1 and a negative effective length is ob- 
tained) . In addition to these changes, it was found that the exponents 
on pressure and final gas velocity obtained for heptane in references 
1 and 2 did not fit all propellants; therefore, new exponents were de- 
termined for these parameters. 

The spread in the effective length required to achieve a given 
percent vaporized was greatest in the 20- to 30-percent vaporized region 
as was also found in references 1 and 2. The variation in the position 
of the inflection point in the mass vaporized curve produced most of 
the spread. The plot of percent -mass unvaporized as a function of ef- 
fective length gave almost a single curve for each propellant. 

Correlated results for various standard deviations of the five 
propellants are shown in figure 3. With all propellants, increasing 
the standard deviation of the spray increased the effective length re- 
quired to vaporize a given high percentage of the propellant. The lines 
in figure 3 represent the log mean effective length for various operating 
conditions. The results of cold-flow spray studies (ref. 5) indicate 
that rocket-engine sprays have a standard deviation of about 2.3; there- 
fore, the results shown in figure 2 should be applicable to most rocket- 
engine systems. This was verified in the comparison of experimental 
and calculated results of reference 2. 
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CONCLUDING REMARKS 

Calculations were made to determine the vaporization rates of liquid 
drops for n -heptane, ammonia, hydrazine, oxygen, and fluorine. The cal- 
culations were made for each propellant with various spray conditions, 
and various engine -design and operating parameters to show how these 
variables would affect the vaporization rate and chamber length required 
to vaporize the drops. The results are correlated by an effective length 
for each propellant. The calculations have shown that the effective 
length required to vaporize a given high percentage of propellant is 
the shortest with oxygen and increases for ^fluorine, n-heptane, ammonia, 
and hydrazine, in that order. 


Lewis Plight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, June 2, 1958 
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APPENDIX A 
SYMBOLS 

surface area of drop,, sq in. 

100 

constant for mass distribution, — - — 

v2jt In cr_ 

e 

coefficient of drag for spheres, dimensionless 
specific heat at constant pressure, Btu/ (lb) (°F) 

diffusion coefficient, sq in. /sec 
heat-transfer coefficient, Btu/ (sq in. ) (sec) (°F) 
coefficient of mass transfer, sec - -*- 
thermal conductivity, Btu/(in. ) (sec) (°F) 
chamber length, in.. 

molecular weight of propellant, Ib/mole 
mass median drop radius, in. 

mass of "i" drop, lb 

mass of "i" drop at beginning of tine, lb 

number of drops in group of "i” sized drops 
chamber pressure, lb/s q in. 
vapor pressure of liquid, Ib/sq in. 
heat-transfer rate of drop, Btu/sec 

percent of mass in drops smaller than "r” 
drop radius, in. 
temperature of gas, °R 
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T c critical temperature of propellant, °R 

T-j temperature of liquid drop, °R 

T r reduced temperature, Tj^ q /T c 

T mean gas temperature, °R 

U velocity difference 'between gas and drop, in. /sec 

u velocity of gas in chamber, in. /sec 

u fin final velocity of gas, in. /sec 

S tfj 

v droplet velocity, in. /sec 

w vaporization rate of fuel, lb/sec 

Z correction factor for heat transfer, dimensionless 

a correction factor for mass transfer, dimensionless 

9 time, sec 

A latent heat of vaporization, Btu/lb 

(i viscosity, lb/(in. ) (sec) 

p density, Xb/cu in. 

0 g geometric standard deviation, ^ 

Subscripts: 
l liquid fuel 

mx vapor mixture 

o initial condition 

p combustion products 

v vapor mixture 


4720 



4720 


NACA TN 3883 


9 


APPENDIX B 
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PHYSICAL PROPERTIES 
Heptane with Gaseous Oxygen 

Density of liquid, 

p z = 3.L662X10" 2 - 9.5355X10“ 6 T z - 6 . 945X10" 9 T 2 , Ib/cu in. 
Heat of vaporization, 

A » 139.9 + 0.181 Tj - 2. 7875X10 " 4 Tj, Btu/lb 

Specific Heat of liquid, 

= 0.231 + 5.62X10“ 4 T 2 , Btu/(lb) (°F) 

Specific heat of heptane vapor, 

c p ^ v = 0.5755 + 1.805Xl0 -4 f, Btu/(lb) (°P) 

Specific heat of combustion products, 

Cp^ p = 0.2898 -f- 4.07X10" 5 T, Btu/(lh) (°F) 

Thermal conductivity of heptane vapor, 

fcy. - 2.914X10 -8 + 5 . 847X10 -11 T, Btu/(in. ) (sec) (°F) 

Thermal conductivity of combustion products, 

kp = 1.3349X10" 7 + 3 . 4111X10 “ 10 T, Btu/(in.) (sec) (°P) 

Vapor pressure, 

lnp = 11.94763 - 5255.89687/^ - 101.58), lb/sq in. 

Viscosity of heptane vapor, 

VLy = 2. 106x10“ 7 + Y^OXlO -10 !, lb/(in.) (sec) 
Viscosity of combustion piroducts, 

Pp = 5.5615X10” 7 + 1.4214X10" 9 f, lh/(in. ) (sec) 
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Diffusion coefficient of vapor mixture, 

D « [-9.815X10' 4 + 1.973X10" 6 T + I.ISIS^O" 9 ^) 2 ]^, sq in. /sec 

Molecular weight -of combustion products,, 

Mp = 31, lb/mole 

Oxygen with Gaseous Heptane 

Density of liquid, 

p z = 0.023079 + 2.7359X10" 4 Tj - 9 .9465x10" 7 t|, lb/cu in. 

Heat of vaporization, 

A = 61.332 + 0.5916 ^ - 2.48xl0" 3 T 9 , Btu/lb 

Specific heat of liquid, 

c^ = 0.3726 + S^SSKLO -4 ^, Btu/(lb) (°F) 

Specific heat of oxygen vapor, . 

Cp v a 0.21333 + 2.2111XL0" 5 T, Btu/(lb) (°F) 

Specific heat of combustion products, • 

c_ — 0.2898 + 4.07Xl0‘ 5 T, Btu/(lb) ( °F) 

Thermal conductivity of oxygen vapor, 

Ky = 2.66UXL0” 7 + 3.4057X10" 10 T, Btu/(in.) (sec) (°F) 
Thermal conductivity of combustion products, 

kp a 1.3349X10" 7 + 3.41HX10" 10 T, Btu/(in. ) (sec) (°F) 

Vapor pressure, 

in p = 11.9584 - H>/b 4 in- 
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Viscosity of oxygen vapor, 

liy. = 2. 2500X10" 7 + 1.1702X10" 9 !, Ib/(in. ) (sec) 

Viscosity of combustion products, 

Up a 5.5615X10“ 7 + 1 . 4214x10 “®T, lb/(in. ) (sec) 

Diffusion coefficient of vapor mixture, 

D = [-2.36396X10 -3 + 5.7897XL0 -6 T + 2. 87685X10 " 9 (T) 2 ] sq in./sec 

Molecular weight of combustion products, 

Mp = 31, lb/mole 


Ammonia with Gaseous Oxygen 


Density of liquid, 

P z =» 3.13955X10" 2 - 6.0219XL0~ 6 Tj - 2. 23169x10 “ 8 t|, lb/cu in. 
Heat of vaporization, 

A a 676.362 + 0.32247 Tj - 1.21Qxl0" 3 Tf, Btu/lh 
Specific heat of liquid, 

e p l a 0.8314 + 5.7993xl0 _4 T2, Btu/(lb)(°F) 

Specific heat of ammonia vapor, 

c p ^ v = 0.60931 + 8 . 3136 1x10 " 8 T, Btu/(lb) (°F) 

Specific heat of combustion products, 

c p ^ p a 0.403578 + 6. 9 96x10 " 6 T, Btu/(lb) (°F) 

Thermal conductivity of ammonia vapor, 

ky. a 4.6853X10" 7 + 3. 551X10 “ 10 T, Btu/(in.) (sec) (°F) 
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Thermal conductivity of combustion products, 

kp = 5. 10692x10”® + 7.616145 x10“ :l 0 T, Btu/(in.) (sec) (°F) 

Vapor pressure. 

In p = 13.1485 - 5 8 . 2351 ) ' in * 

Viscosity of ammonia vapor, 

= 1.07097x10”® + 8.11607xl0” 10 T, Tb/(in. ) (sec) 

Viscosity of combustion products, 

= 1.64352x10" 7 + 2. 45 25 2x10” 9 T, lb/(in. ) (sec) 

Diffusion coefficient, 

D - (-4. 87304X10“ 2 + 4. 01198x10 - 8 .358X10“ 10 T 2 ) sq in. /sec 

Molecular weight of combustion products,- 

Mp = 24, lb/mole 


Hydrazine with Gaseous Oxygen 

Density of liquid, 

P z = 3.06 2318x10” 2 + 4.028897X10" 5 T Z - 5 .54321 x10“ 8 t|, lb/cu in. 
Heat of vaporization, 

A =» 730.747 - 0.3591305 T^ + 1. 214XL0” 4 T^, Btu/lb 
Specific heat of liquid, 

c p,2 ■ 0.589125 + 2.80708Xl0“ 4 Tj, Btu/(lb) (°F) 

Specific heat of hydrazine vapor, 

c pjV = 0.3360 + 1.804X10" 4 T, Btu/(lb)(°B) 
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Specific heat of combustion products , 

Cp^ p = 0.403578 + 6.996xlO“ 6 T, Btu/(lb) (°F) 

Thermal conductivity of hydrazine vapor, 

kV = 1.23753X10" 8 + 2.23O358xl0- io T, Btu/(in. ) (sec) (°F) 

Thermal conductivity of combustion products, 

kp = 5.10692X10" 8 + 7 . 616145xlO -10 T, Btu/(in. ) (sec) (°F) 

Vapor pressure, 

1“ p - 14 - 328787 - ■ ( Il 73 - 6 6^17X5 ;'> to- 

Viscosity of hydrazine vapor, 

JV = 4.19981XL0" 8 + 9.581164X10- 10 f, lb/(in.) (sec) 

Viscosity of combustion products, 

Mp = 1.6435 2x10"" ^ + 2.45252X10“ 9 T, lb/(in. ) (sec) 

Diffusion coefficient of vapor mixture, 

D = (-5.3537X10" 4 + 3.13874XL0 -6 T + 3. 37045x10 _9 T 2 ) , sq in. /sec 

Molecular weight of combustion products, 

Mp = 24, Ib/mole 

Fluorine with Gaseous Hydrogen 

Density of liquid, 

Pj = 6 . 846x10“ 2 - 3.036xl0" 5 Tj - 3.9308 x10" 7 t|, Ib/cu in. 
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Specific heat of liquid, 

c p ^ = 0.349 + 1.212xlO“ 4 T 1 , Btu/(lb) (°F) 

Specific heat of fluorine vapor, 

c = 0.223994 + 1.667XlO” 6 f, Btu/(lb) (°F) 

t? J V 

Specific heat of combustion products, 

c » 0.32332 + 2. 06 8x10 "®T, Btu/(lb) (°F) 

1 ? ) 3 ? 

Thermal conductivity of fluorine vapor, 

\ - 2.7606X10' 7 + 3.47122xlO“ 10 T, Btu/(in.) (sec) (°F) 

Thermal conductivity of combustion products, 

kp = 1.0765X10 -7 + 4 . 31280x10 “ 10 T , Btu/(in.) (sec) (°F) 


Vapor pressure, 

In p a 12.3171 


1482.9545 



Tb/sq in. 


Viscosity of fluorine vapor, 

Py =« 1.2591x10”® + 1.584X10" 9 !, lb/(in. ) (sec) 

Viscosity of combustion products, 

Pp = 2. 6009 8x10" 7 + 1 . 0417 8X10”® T , Ib/(in.) (sec) 

Diffusion coefficient of vapor mixture, 

D = [-4.5732X10" 3 + -1.078x10 _5 T + 5.421X10" 9 (f) 2 ] sq in. /sec 

Molecular weight of combustion products, 

Mp = 20, lb/mole 
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Figure 1. - Continued, typical droplet histories for various propellant a (see table X for conditions). 
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Figure 2. - Correlated results for propellant with standard deviation of 2.3 
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Figure 2. - Continued. Correlated results for propellant with standard deviation of 2.3. 
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Figure 2. - Continued. Correlated results for propellant vith standard deviation of 2.3. 
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(i) Mass vaporized for liquid-fluorine sprays. 

Figure 2. - Continued. Correlated results for propellant with standard deviation of 2.3. 
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Figure 2. - Conc lude d - Correlated results for propellant with standard deviation of 2.3. 
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